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Abstract. Most of the Skyrme interactions are known to predict spin or isospin 
instabilities beyond the saturation density of nuclear matter which contradict 
predictions based on realistic interactions. A modification of the standard Skyrme 
interaction is proposed so that the ferromagnetic instability is removed. The new 
terms are density dependent and modify only the spin p-h interaction in the case of 
spin-saturated system. Consequences for the nuclear response function and neutrino 
mean free path are shown. The overall effect of the RPA correlations makes dense 
matter more transparent for neutrino propagation by a factor of 2 to 10 depending of 
the density. 
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The understanding of dense nuclear matter and its application to nuclear physics 
and to astrophysics is strongly related to the properties of in-medium nuclear interaction. 
During the recent years, much efforts went into the fundamental understanding of the 
bare nuclear interaction, either through the developement of effective field theories 
derived from QCD symmetries [I] or through renormalisation group arguments applied 
to the Lippmann-Schwinger equation which generate a phase-equivalent low energy bare 
interaction named t>i ow -fc [2J. These theories aim at a separation of the long-range 
component of the nuclear forces dominated by the pion exchange, which is well under 
control, from the intermediate and short range components, which are dominated by 
correlated pion and heavy meson exchanges that are poorly known. A close agreement 
between the interaction tw-fc and the effective Gogny interaction [3j have been found 
concerning the predictions of pairing properties in nuclear matter in a wide range of 
sub- nuclear densities [4]. 

Another modelisation of the nuclear interaction in dense matter based on the 
symmetries of QCD, the quark-meson coupling model, see Ref. [5] and references therein, 
has given a sound foundation to the Skyrme effective nuclear interaction [6l Ej. The 
modern Skyrme forces are commonly adjusted to the empirical knowledge of nuclear 
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matter at saturation density, and sometimes include also theoretical calculations of 
asymmetric or pure neutron matter, as well as a set of nuclei described within Hartree- 
Fock + BCS or Hartree-Fock Bogoliubov theories. So far the Skyrme interactions are 
successfully used to describe properties of nuclei such as binding energies, radii and 
excited states. Nevertheless, due to its very simple functional form, the standard Skyrme 
interaction is not suited for being extrapolated in astrophysical situation such as neutron 
stars. Indeed, the mean field solution obtained with such interactions become unstable 
for densities larger than the saturation density po and also for neutron rich matter. 

A recent extensive study of the symmetry energy deduced from Skyrme 
interactions [E] has shown that the symmetry energy becomes negative beyond the 
saturation density of nuclear matter for many of these interactions. This is related to 
the instability of the mean field with respect to the isospin density fluctuation Spt where 
Pt = Pn — Pv Depending on the parameters of the interaction, instabilities can occurs in 
different channels, at various densities and at different isospins. It has been shown that 
there is a limited interval of densities between po and 3po for which the mean field is 
stable in symmetric matter and also in neutron matter The stability of the ground 
state toward small fluctuations can be studied within the HF+RPA framework. The 
fluctuations around the mean field are induced by the particle-hole interaction, 



in terms of the dimensionless Landau parameters F#, Fj,, G#, G' e and Legendre 
Polynomials Pi(cos8). N = gm*kp /(27r 2 'H 2 ) is the density of state around the Fermi 
energy and g is the degeneracy. The matter is stable unless one of these parameters of 
multipolarity I becomes lower than —21 — 1. 

Recently, the ferromagnetic instability has regained interest [TOj HI] due to the 
observation of extremely high magnetic fields (~ 10 15 ~ 16 G) in compact stars [12] and 
also of giant flares observed recently on 27 December 2004 [13]. We show in Fig. [I] the 
critical density pf at which the asymmetric matter becomes unstable with respect to 
the spin fluctuations. The proton fraction x p = p v j p is changed from symmetric matter 
(x p = 1/2) to neutron matter (x p = 0). We select various Skyrme interactions which 
are commonly used in the description of finite nuclei. Depending on the interaction, the 
density p f can be very close to p (BSkl6 [14], RATP [15]), around 2p (SkM* [IB] . 
SLy5 [I?]) or nearly reach 3po (SGII [IE], LNS [19]) in symmetric nuclear matter 
(x p =l/2). When the proton fraction is decreasing, the density p f either decrease (BSkl6, 
RATP, SkM*, SGII) or increase (SLy5, LNS). From Fig. [Hit could be deduced that the 
prediction of pf is varying largely among the selected Skyrme interactions. 

According to microscopic calculations with realistic interactions such as diffusion 
Monte-Carlo [20] or Brueckner Hartree-Fock (HF) calculations [2Tj [22| [23] . there is 
no ferromagnetic instability up to substantially high densities. The onset of the 
ferromagnetic instability beyond a few times po represents a serious limitation of the 
Skyrme interaction. This limitation should be cured in the calculation of nuclear matter 
properties such as equation of state, the response function and microscopic processes 
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such as neutrino mean free path. It is then timely to enrich the effective Skyrme 
interaction in order to extend its domain of application and to take advantage of its 
simple form. To this end, extended Skyrme interactions shall mimic more accurately 
microscopic results of realistic interactions and preserve the accuracy of the original 
Skyrme interaction in the description of nuclei. 

It is known that the spin and the spin-isospin excitations have strong impact on 
astronomical observables such as neutrino mean free path in neutron star, Oz/— and 2v 
double beta decay processes. On the other hand, it has been recognized also that Skyrme 
interactions have a serious shortcoming in the spin channels. Our main purpose in this 
manuscript is to improve the spin dependent parts of Skyrme interactions keeping its 
simplicity and good properties for ground state properties. In this way, we will be able 
to extend the field of possible applications of Skyrme interactions for spin dependent 
excitations not only in finite nuclei but also infinite nuclei within the self-consistent 
theoretical model. It is important to make a bridge between finite and infinite systems 
by the self-consistent model without introducing any extra parameters. Top of realistic 
spin and spin-isospin interactions, we would like to introduce further tensor interactions 
and two-body spin-orbit interactions for constructing a global energy density functions 
of spin and spin-isospin channels. In this article, we first propose a simple way to get 
ride of the ferromagnetic instability by introducing a limited number of new terms to 
the standard Skyrme interaction in Sec. [U In Sec. [2} we apply the new interaction to 
the calculation of the RPA response function in nuclear matter and then we calculate 
the neutrino mean free path. Conclusions are given in Sec. [3j 

1. Extended Skyrme interaction 

Some of the most recent Skyrme interactions are fitted so as to reproduce the theoretical 
binding energy in symmetric and neutron matter up to a density as much as few times 
Po. As a result, the instability in the isospin channel, which occurs only if the binding 
energy of neutron matter becomes more attractive than that of symmetric matter, is 
usually either totally removed or pushed to very large densities. It is then possible to get 
ride of the isospin instability without modifying the standard Skyrme interaction. In the 
following, we focus on the spin instability and test the new terms added to the existing 
Skyrme interactions which do not show the instability in the isospin channel. We adopt 
three different Skyrme interactions, SLy5, LNS and BSkl6 for this study. SLy5 is an 
interaction suited by construction to make predictions in neutron rich nuclei, LNS is an 
interaction which reproduces global features of G-matrix in symmetric and asymmetric 
nuclear matter and BSkl6 reproduces known masses of nuclei with the best accuracy of 
rms deviation 632 keV. Some nuclear matter properties of these interactions are listed 
in Tab. HJ The monopole Landau parameters in symmetric (SM) and neutron (NM) 
matter are drawn in Fig. [2] as a function of the density. As expected, these interactions 
has no instability in the isospin channel. Notice however that in the isospin channel, 
the Landau parameters Fq sm obtained from the three interaction are more attractive 
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than the one deduced from the microscopic G-matrix calculation [23] (filled squares in 
Fig. [2]), especially in the high density region. However, it is not necessary to correct this 
channel since there are no instabilities in the density range considered in Fig. [2j In the 
spin channels, the Landau parameters Go,sm, ^?osm an< ^ ^o,nm become more attractive 
as the density increases. The three Skyrme interactions have pathological behaviors in 
the spin channels G and G' Q for which we now propose a prescription to cure. 

From the density dependence of the Landau parameters represented in Fig. [2j it is 
clear that the standard Skyrme interactions are not repulsive enough in the spin channel, 
and most probably also in the isospin channel at higher density than the saturation 
density. The repulsion shall come from the effect of the three body force [21]. Through 
the density matrix expansion, it has been shown that in-medium many-body correlations 
give rize to density dependent terms in the nuclear functional [25J . These terms could be 
directly included in an effective density dependent nuclear interaction such as Skyrme 
or Gogny. The most important contribution comes from the scalar-density and it is 
the only density which has been considered. However, spin, isospin and spin-isospin 
density dependent terms shall also be considered from the density matrix expansion. 
Here, we will study the impact of the new density dependent terms such as the spin and 
spin-isospin densities: 

K add -( ri ,r 2 ) = ±t s 3 (l + x s 3 P a )[p s (R)r°8(r) + +xi t P a )[p st (R)]'r-8(v) 

(2) 

where P a = (l+cri-0"2)/2 is the spin-exchanged operator, r = ri — r 2 and R = (ri+r-2)/2. 
In Eq. ([2]), we have introduced the spin density p s = P]— p\ and the spin-isospin density 
Pst = Pni — Pni — Ppt + Pvi- Spin symmetry is satisfied if the power of the density 
dependent terms 7 S and j st are even. 

In the literature, it has already been proposed to add new terms to the standard 
Skyrme effective interaction. One of the motivation is that the HF single particle levels 
are usually too widely spaced at the Fermi surface. This shall be corrected by including 
the coupling of single particles to vibrations which increase the effective mass of the 
single particle states from m*/m < 1 to m*/m > 1 at the Fermi surface. This dynamical 
effect could be mimic by introducing a gradient term in the density dependent interaction 
(see Refs. [26], [27] and references therein). Such a term has also been included in global 
fits of mass formula to make nucleosynthesis calculations as accurate as possible [28J. 
However, the removal of the spin instability makes symmetric matter unstable at high 
density [29J. Up to now, there is no satisfactory additional term which prevent the 
matter to fall into spin instabilities. 

In the following, we adopt the notations of Ref. [T7J where the density functional TC 
is expressed as a sum of the kinetic term /C, a zero range term Tio, a density dependent 
term 7i 3 , an effective-mass term 7i e fj and some additional terms coming from spin-orbit 
coupling, spin and gradient coupling and coulomb interaction. The additional terms (T2]) 
modify the density dependent part, 7i 3 to be T^+T^+T^*, where the additional density 
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nf = f 2 p:t [(i + fV + x fpl - (*? + + p 2 ,) - ~(pL + A)] (4) 

where p sn = p n f — p n ± and p sp = p P i — p p ±. The mean field U q , where q — n,p, is then 
corrected to be 

Uf A - = §P?{(2 + xl)p - (1 + 2xl)p q } + {(2 + *f )p - (1 + 2*f )pj. 

(5) 

In symmetric nuclear matter the Landau parameters Go and G' [30] are also modified 
by the following additional terms 

/^add. -is -ist i 

V = l 7 ^ 7 * - Wpr 2 - (2^ + - pi] + g(xf - \) P :t 



+ |(x3-^)(7, + l)(7 s + 2)p?, (6) 

/'"Vadd. j.s 

= ^(%t - l)[3p 2 + (2x^ - 1) p\ - (2s? + l)pMr 2 - f A Pl s 

-%hst+2)( lst +i) P :r . (?) 

In pure neutron matter, the additional terms for Go are 

/'"fadd. js j.st 

^- - % - 4)7,(7, - l)pVr 2 + ttt (af - 1) P 7st 



iV n 24 v J/ ' v ' " 1 s 12 



+ f (*S-l)(7. + 2)(7. + l)p?. (8) 

Non-trivial result for the Landau parameters Gg dd ' and G ad<1 imposes a condition 
on the powers of the density dependence in Eq. (j2|) to be uniquely fixed as 7 S = 2 and 
7 st = 2. Notice that spin symmetry is also satisfied by these selections. It turns out 
that the additional contributions (jSJ) to the mean field are simply null in spin-saturated 
systems, as the densities p s = and p st = 0. It is thus possible to add the new terms (j2D 
to the existing Skyrme forces without destroying good properties of the original Skyrme 
interactions. We obtain the contributions to the Landau parameters 

or = ft^ , o) 

G' ad± = f^P 2 , (10) 
in spin-saturated symmetric matter and 

Gf d - = ^m-xt) P 2 . (ii) 

in spin-saturated neutron matter. The contributions to the Landau parameters Go and 
G' in symmetric matter in Eq. (Q and (jTUl) depend only on the parameters i| and tf, 
while that for the Landau parameter Go in neutron matter in Eq. ( fTTT) depends on t| 
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and X3. We shall then first fix i| and tf so as to reproduce the Landau parameters in 
symmetric matter, and then we fix X3 in neutron matter. Since the Landau parameters 
Go an d G' are independent of xf, this parameter can be set as xf=0. 

We show in Fig. [3] the contributions of the new terms for the Landau parameters Go 
and G in symmetric matter and for Go in neutron matter added to the three effective 
interactions SLy5, LNS and BSkl6. We compare the original Landau parameters in 
Fig. [2] (solid lines) with the ones which include the contributions of the new terms ([6]), 
(J7]) and (jSJ). The dashed lines correspond to different values of the new parameters 
£3, tf and X3 with the step indicated in each graph. The values for the parameters 
i|, tf and x% resulting from the best adjustment for asymptotic behavior at very high 
density to the Brueckner HF calculations are given in Tab. [2] and are labelled SLy5st, 
LNSst and BSkl6st respect ivelly. With the parameters given in Tab. [21 it is confirmed 
that the correction terms to the Landau parameters (j91)- (fTTl) are repulsive enough to 
remove the spin instabilities. There are however still differences between the corrected 
Landau parameters and the ones given by the G-matrix which could be originated from 
the tensor interaction. We can see that the LNS interaction, fitted originally to the 
equation of state deduced from the G-matrix, gives the best fitted Landau parameters 
among the three interactions. Finally, it could be noticed that the corrected Landau 
parameters G'q™ are increased by about +0.3 which makes it closer to the empirical 
value. Notice that the empirical value is estimated to be 1.3±0.2 from Wood-Saxon 
single-particle states plus one-pion and rho meson models (see for instance Tab. I of 
Ref. [3T] and references therein and also Ref. [32]) and 0.7±0.1 from self-consistent 
Skyrme mean- field models [33J. 

2. RPA response functions and neutrino propagation 

As a consequence of the new terms the spin channels are more repulsive than the 
original Skyrme functional. The effect of the new terms could then modify the response 
function of collective spin modes already at po- For instance, we show in Tab. [2] the 
Landau parameters G and G deduced from the original Skyrme interactions and the 
one including the new terms G^ ew and G'^ ew . The increased of the Landau parameters 
may have significant influence on spin and spin-isospin excitations such as Gamow- Teller 
states. 

In the following, we calculate the RPA response function in nuclear matter at finite 
temperature, 

5^(g, W| T) = -- , * W'^W.T) (12) 
7r 1 — e ^i 1 

where % '^'(ftw,^ is the susceptibility [31] obtained as the solution of the Bethe- 
Salpeter equation [35]. In Fig. (Ill) , we show the response functions in the spin channels 
(S = 1) at T = MeV and at densities p = po and p = 2p calculated by the LNS 
and LNSst interactions. The HF solution (dotted line) is compared with the RPA using 
the original LNS Skyrme interaction (dashed line) and also with the RPA including 
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the new terms in LNSst (solid line). According to the semi-classical Steinwedel- Jensen 
model [36] the optimal transfered momentum to compare with nuclei shall be q = tc/2R 
where R is the radius of the nuclei. For 208 Pb, one thus obtain g=0.22 fm _1 . For p = po, 
the effect of the new terms is to move the collective mode to slightly higher energies by 
0.5-1 MeV. For p = 2po, the (S — 1, T — 1) Gamow- Teller channel is not far from being 
unstable from the original Skyrme interaction and the low energy collective mode is 
being formed. By the new term tf, the low energy mode is suppressed and the strength 
is reduced substantially at u > 0. A reduction of the strength at low energy is also 
observed for the (S = 1,T = 0) channel. This effect shall be seen in the calculation of 
microscopic processes such as the neutrino cross section. 

Neutrinos play a crucial role in physics of supernova explosions [37] and in the 
early evolution of their compact stellar remnants [381 [39j HQ]. To clarify the effect 
coming from the additional interaction (T5]), we consider the case of pure neutron matter. 
The scattering of neutrinos on neutrons is then mediated by the neutral current of the 
electroweak interaction. In the non-relativistic limit, the mean free path of neutrino 
with initial energy E u is given by [4"T] 



where Gf is the Fermi constant, cy (ca) is the vector (axial) coupling constant, 
k\ = (E u , ki) and k% are the initial and final neutrino four-momenta, while q = k\ — k% 
is the transferred four-momentum, and cos# = ki • k3. In the following, we impose 
the average energy E v = 3T in MeV [40J. In Eq. ffl3|) . the contribution of the 
response function S^ s \q,u,T) is clearly identified. It describes the response of neutron 
matter to excitations induced by neutrinos, and contain the relevant information on 
the medium. The vector (axial) part of the neutral current gives rise to density (spin- 
density) fluctuations, corresponding to the spin S — (S — 1) channel. 

The neutrino mean free path is shown in Fig. as a function of the baryonic 
density and for different temperatures (T=l, 10, 20 MeV). The neutrino mean free path 
calculated with the original LNS Skyrme interaction (dashed line) is strongly reduced at 
high density as a consequence of the onset of the spin instabilities |4"3~1 142] . The inclusion 
of the new density dependent interaction (T5]) in LNSst removes the spin instability and 
reduces the response function as shown in Fig. HI As a consequence, the mean free path is 
increased and the matter is more transparent to neutrino as the density increases. This 
result is similar to a previous calculation of neutrino mean free path deduced directly 
from a microscopic G-matrix |44) . It illustrates the important contribution of the spin 
channel to the neutrino mean free path at high density. It is also clear from Fig. [5] that 
the effects of the correlations are not washed out by the temperature in the range going 
up to T = 20 MeV. 




(13) 
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3. Conclusions 



We have proposed the extension of the Skyrme interaction (j2J) which conserves its 
simplicity and extends its domain of stability to large densities and large isospin 
asymmetries. The parameters of the new terms have been adjusted based on microscopic 
G-matrix calculations. For the three interactions SLy5st, LNSst and BSkl6st, the spin 
channels become more repulsive compared to the original interactions SLy5, LNS and 
BSkl6; and the dimensionless Landau parameter G' is increased by about +0.3. The 
simplicity of the extended Skyrme interaction makes possible extensive calculations not 
only in asymmetric and dense matter but also in asymmetric finite nuclei. We have 
calculated RPA response functions of spin channels and also neutrino mean free path in 
dense matter. The collapse of the neutrino mean free path at high density is suppressed 
by the newly added terms in the Skyrme interaction and the overall effect of the RPA 
correlations makes dense matter more transparent for neutrino propagation by a factor 
of 2 to 10 depending on the density. 

We would like to extend further the applications of our new spin-density dependent 
interactions for the systematic study of odd mass nuclei, the spin-dependent excitations 
in finite nuclei and also infinite nuclear matter. Eventually our ultimate goal is to 
construct the global energy density functional for the spin and spin-isospin dependent 
interactions including tensor and two-body spin-orbit forces. Concomitantly, an 
extensive experimental campaign dedicated to the measurement of Gamow- Teller and 
also other spin-isospin modes is proposed at RIKEN [35]. The new data will be very 
usefull to compare with our global energy density functional. 
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Table 1. Nuclear matter properties of the Skyrme interactions SLy5, LNS, BSkl6 at 
saturation density po; the binding energy a v , the incompressibility K^, the symmetry 
energy a s , the density pf of the ferromagnetic instability. 





Pa 


a v 




a s 


Pf/Po 




(mi" 3 ) 


(MeV) 


(MeV) 


(MeV ) 




SLy5 


0.16 


-15.97 


229.9 


32 


2.09 


LNS 


0.1746 


-15.32 


210.9 


33.4 


2.48 


BSkl6 


0.1586 


-16.05 


241.6 


30 


1.24 



Table 2. Parameters of the additional terms i| (in MeV.fm 37s ~ 2 ), tf (in 
MeV.fm 37st_2 ) and for the interactions SLy5st, LNSst and BSkl6st in order to 
reproduce the realistic Landau parameters by Brucckncr HF calculations. The powers 
of the density dependence in Eq. (2) are set to be j s = j st — 2, and x|* = 0. Arc 
also shown the dimensionless Landau parameters Go, Gq cw , G' and G' new deduced 
from the original Skyrme interactions and the new interactions including the additional 
terms, respectively. 
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Figure 1. Ferromagnetic phase diagram for various effective Skyrme interactions. The 
horizontal axis shows the critical density pf at which the asymmetric matter becomes 
unstable, while the vertical axis shows the proton fraction x p = p p /p. The matter is 
spin symmetric for smaller density than pf, while it becomes ferromagnetic for larger 
density than pf. See the text for details. 
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Figure 2. Landau parameters in symmetric nuclear matter (SM) and neutron matter 
(NM) as a function of the density p obtained by the Skyrme interactions and the 
Brueckner HF calculations using 2BF+3BF Solid line for SLy5, dotted line for 

LNS, dashed line for BSkl6 and filled squares for Brueckner HF calculations. 
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Figure 3. Landau parameters for spin and spin-isospin channels Go and G' in 
symmetric nuclear matter and neutron matter. The solid curves are the original ones, 
while the filled squares are obtained by the Brueckner HF calculations as shown in 
Fig. [21 The dashed curves correspond to different values of the parameters i§, if* and 
x§. The parameters are changed from bottom to top multiplying integers (f ,2,3,-) by 
the mesh size given in each window. The Landau parameter Go,nm is calculated with 
the optimal value for i§ in Tab. El which reproduces best the Brueckner HF Go,sm 
values. The thin dotted lines in the bottom panels correspond to the results with 

4=0. 
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Figure 4. RPA response functions S^ S ' T \q,u!) for q=0.22 fm^ 1 with the temperature 
T=0 MeV calculated with LNS and LNSst interactions. Top panels are calculated for 
(S,T) = (1,0) and (1, 1) with p = po and bottom panels are for the same channels 
with p — 2po. 
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Figure 5. Neutrino mean free path in neutron matter for different temperatures 
T=l, 10 and 20 MeV. The neutrino energy is set to be E v = 3T. The mean free path 
is calculated with the LNS Skyrme interaction by the HF mean field approximation 
(dotted line), including the RPA correlations with the original LNS interaction (dashed 
line) and RPA with the modified spin channel in LNSst (thick line). 



